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ABSTRACT 



In this note we present complete, closed-form expressions for random relative velocities between colliding particles of arbitrary size in 
nebula turbulence. These results are exact for very small particles (those with stopping times much shorter than the large eddy overturn 
time) and are also surprisingly accurate in complete generality (that is, also apply for particles with stopping times comparable to, or 
much longer than, the large eddy overturn time). We note that some previous studies may have adopted previous simple expressions, 
which we find to be in error regarding the size dependence in the large particle regime. 
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1. Introduction and outline 

Gas in astrophysical environments is often in a turbulent state of 
motion, constantly affected by temporally and spatially varying 
accelerations from eddies having a variety of scales. A particle, 
due to its inertia, does not instantaneously follow the gas mo- 
tions but requires a certain time in order to align with the gas 
motion. The particle's interaction with the gas is captured in the 
definition of the stopping time of the particle (sometimes also 
referred to as friction time), 



t, = 



4cgp g o-' 



(1) 



where c g and p g are, respectively, the sound speed and the vol- 
ume mass density of the gas, and m and <x the mass and pro- 
jected surface area of the particle. Due to this inertial lag, a par- 
ticle develops a relative velocity with respect to the gas. In addi- 
tion, these lags also cause particles to acquire relative velocities 
among themselves. 

While the general problem of calculating these relative ve- 
locities has received consid erable attention i n the basic fluid 
dynamics community (see ICuzzi & HoganI d2003l) for refer- 
ences; henceforth CH03), the formalism most frequ ently used 
in the astro physics communit y was developed by IVolk et al.l 
d!980h and iMarkiewicz et all d!991l) (henceforth MMV). In 
these works the final results are given in terms of integrals that 
were not solved analytically. Some workers have used simple 
fits to these numerical results in their models of dust coag- 
ulation; however, simple closed-form expressions for particle- 
particle relative velocities wou l d help streamline these mod- 
els (e.g. ISuttner & Yorkel 120011: lDuTlemond& Dominikl 120051; 
iNomura & Nakagawd2006tlOrmel et al.l2007l) . Recently, CH03 
obtained closed-form expressions from the MMV model for par- 
ticle velocities in inertial space V p , for particle-gas relative ve- 



locities V pg , and for relative velocities between two identical par- 
ticles V pp , but did not extend their results to the general case 
of two particles of different stopping times. Moreover, CH03 
stressed the validity of their analytical results for particles with 
stopping times much shorter than the large eddy turnover time. 
In this note we generalize the approach and results of CH03 to 
obtain closed-form expressions for relative velocities between 
particles of arbitrary, and unequal, size. In Sect. [2] we define 
important quantities and review previous work. In Sect. [3] we 
present two independent approaches for obtaining the desired 
closed-form solutions. In Sect. [4] we give our conclusions and a 
summary. 

2. Definitions and previous work 

Nebula gas turbulence is generally described as being composed 
of eddies having a range of spatial scales I and spatial frequen- 
cies k = l/£, with an energy spectrum E(k) cc k~ 5 ^ and total 
energy Vg/2 per unit mass providing the normalization condi- 
tion 



dk E(k), 



(2) 



from which E(k) = Vg/3&L ■ The largest, or integral 

scale, eddies have spatial scale L = 1 jki_, and the smallest, or 
Kolmogorov scale, eddies have spatial scale r\ = l/k v . The form 
of E(k) given above is the inertial range exp ression most often 
assumed, with E(k) = for k > k n or k < h . IVolketall 0980) 
used a spectrum P(k) = 2E(k) and stipu l ated n o smallest scale n 
for the turbulence, but lWeiden schilling ( 1984) and MMV noted 
that a finite value for n > had profound effects on the par- 
ticle velocities, especially the relative velocities V pp for small 
particles. Each eddy wavenumber k has a characteristic velocity 
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V(k) = y/2kE(k) and overturn time t k = i/V(k) = {W{k)T l . Our 
standard definition of the particle Stokes number is St = ? s /?l, 
where ?l is the overturn time of the largest eddy, generally taken 
to be the local orbit period. The local turbulent intensity is de- 
scribed by its Reynolds number, Re, defined as the ratio between 
the turbulent and the molecular kinematic viscosities, Re = vt/v. 
The values for (, v and t at the integral scale then follow from 
Re, e.g., r] — Re~ 3 ^ L and t n = Re' 1 ^ 2 ^. These expressions bring 
Re into the final expressions for particle velocities as a limit on 
certain integrals (cf. CH03 for more detail). In the notation of as- 
trophysical "cr-models", Re = ac„H g /v = ac 2 /vQ. where c s , H s , 
and v are the sound speed, vertical scale height, and kinematic 
viscosity of the nebula gas and Q is the orbital frequency. 

IVolk et al.l d 1980h introduced the concept of "eddy classes". 
Class I eddies vary slowly enough that a particle, upon enter- 
ing a class I eddy, will forget its initial motion and align it- 
self to the gas motions of the eddy before the eddy decays 
or the particle leaves the eddy. Class II eddies, on the other 
hand, have fluctuation times shorter than the particle's stop- 
ping time ? s , and fluctuate too rapidly to provide more than a 
small perturbation on the particle. The timescale on which an 
eddy decays is given by tk, while the eddy-crossing timescale is 
fcross ~ £/Vre\ = {kV K \(k)Y l , with V re i the relative velocity be- 
tween a grain and an eddy. For an eddy to be of class I both 
and f cross must be larger than the particle's stopping time. The 
boundary between these classes occurs at k = k* (or at tk = f ) 
which can be defined as dVolketalJ 0980), MMV): 



and the same analytical result is obtained, i.e. CH03, 



1 _ 1 

t, ~ t* t. 



— = -+k t V Te i(k*). 

ross * 



(3) 



It is important to realize that k* (or f* ) is a function of stopping 
time f s , that is, the boundary separating the two classes is dif- 
ferent for each particle. The different treatment for the two eddy 
classes k < k* and k > k* forms the core of the derivation of the 
turbulence-induced particle velocities. 

All turbulent velocities in this note are statistical, root-mean- 
square, averaged quantities. The average inertial space particle 
velocity V p is given by Eq. (6) of MMV. 



dk2E{k)[\ -K 2 ) 



+ f" dk2E(k)a-K)[g(x)+Kh(x)], (4) 

in which K = t s /(t s + tk). The K 2 term in the first integral 
results from the more recently preferred "n = 1" gas veloc- 
ity autocorrelation function (MMV and CH03). The functions 
g(x) = ^tan- 1 ^) a nd h(y) = 1/(1 + y 2 ) with x = Kt k kV K i 
were first obtained bv lVolket ail 0980). 

CH03 noted that, for very small particles with f s <sc ?l or 
St <sc 1 , the second integral becomes negligible, leaving only 
the first integral which is analytically solvable and for which the 
upper limit can be extended to k l; with negligible error. Here, to 
generalize the approach of CH03 to particles of arbitrary size, 
we approximate h(x) = g(x) = 1 f° r a ^ particle sizes (see 
CH03 Sect. 2.2.3 for supporting logic). Numerical calculations 
of h(x) and g(x) validate this approximation to order unity (see 
Appendix[A]i, and we gain further confidence in it from a poste- 
riori comparison with exact numerical model results. The gen- 
eral expression for V? is then the same as in the f s «: ?l regime, 



r-K, 

V 2 = dk2E(k)[\-K 2 ) 

Jk L 

St 2 (l - Re- 1 ' 2 ) 



V? 1 



(St + \){$>l + Re- l l 2 ) 



(5) 
(6) 



CH03 did not give this explicit result for Vp, but merely noted 
that it was straightforward to derive it from their Eq. (19) for V p£ , 
and the general relationship V 2 = V 2 - V 2 ; however we will use 
it explicitly here. 

Comparison of the predictions of this simple expression with 
detailed numerical results (MMV, CH03) show that it is indeed 
a good approximation for arbitrary St. A more accurate approx- 
imation to Eq. (O, in which the g and h functions are approxi- 
mated as power-laws in k*/k, is outlined in AppendixlAl Unless 
St <K 1, we can neglect the Reynolds number term in Eq. (|6]l and 
obtain V D = VJ Vl + St, a well known result ( TVolk et aljfl980l: 
ICuzzi et alj|1993t ISchrapler & Henningll2rj04h which describes 
the diffusivity of large particles in turbulence. 



3. Results 

3.1. k-space approach 

MMV (their Eq. (7)) expressed the relative velocities V p i P 2 be- 
tween particles of different stopping times t\ and t^ as 



V: 



plp2 



2V pI V p2 s AV z l2 . 



(7) 



Having already derived (i — 1, 2) above, we can determine 

AV12 by evaluating the cross term V p i V P 2; this paper presents an- 
alytical solutions of this problem obtained in two separate ways. 
In this subsection we retain the wavenumber dependence; in 
the next subsection we transform to time variables. In Eq. (8) 
of MMV the cross term is given as a sum over the two par- 
ticle sizes involved, which we separate here, writing AV 2 2 = 
V£ + V£-(V* +V«), where 



V: 



-—if 



E(k)dk 

~mi 



r-min(k' r k' 2 ) I 1 \ 

I ' e «(tt«) 



2 \ 

dk 



(8) 



Changing variable to x = k/k^, substituting for E(k), and con- 
verting stopping time f,- to Stokes number St,- = ti/t-c. 



Vl = 



2V ^ 
3(fi + ti) 



j\- 5/3 dx-£ 



St 2 dx 



x 5 / 3 (St,- + x-2/ 3 ) 2 



, (9) 



where we have taken, without loss of generality, k\ < k* v The 
first integral is trivial and the second integral can be solved ex- 
actly as in Eqs. (17-19) of CH03. In evaluating the specific value 
of the integrals above, we need a closed form for the upper limit 
x\ = kyk^. A simple prescription is readily found by inspec- 
tion of Fig. 3 of CH03: x\ = k*/k h = 0.5St^ 3/2 + 1. That is, 
the boundary eddy for particles with stopping time t\ is that for 
which tk ~ t\ until t\ > t^, beyond which it remains constant. 
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This is merely a convenient mathematical shorthand to keep ev- 
erything in closed form. Then, repeating the analytical solution 
of CH03 (Eqs. (17-19)) we obtain 



where the (1 <-> 2) symbol indicates interchange between par- 
ticles 1 and 2. With further manipulation and cancellation of 
terms, the previous expression simplifies slightly to 



Vi = V: 



ti 



E (?i + h) 



[» 



«-2/3n 



St; 



St, 



1 + St, 1 + St;jc* 



2/3 



(10) 



This solution for the cross term is easily combined with 
Eq. © to obtain expressions for particle-particle relative veloci- 
ties AV p. Further manipulation of these expressions may be pos- 
sible, but the important point here is that AVn can be expressed 
in closed form as function of Sti, Sta, V g , and Re. With a few 
minutes of algebra, simpler expressions can be found in the lim- 
iting regimes of interest (Sti <K 1,» 1, etc.) which agree well 
with those which we present in the next section, for analytical so- 
lutions obtained in the time domain instead of the wavenumber 
domain, and where an analytical solution for the boundary k*(f) 
is used rather than the form for x\ adopted above. It should be 
recalled that, for very small particles f s < t n , x* has an upper 
limit of kr,/k L = Re 3/4 (see, e.g. CH03 Fig. 3). 

3.2. t-space approach 

The integrals expressing Vl i and are transformed into a 
simpler form by changing variables from k to t k . Since tk = 
l/kV(k) = [k V2Ar£(/t))~' and E(k) = Ak~ 5/3 for a Kolmogorov 
power spectrum (where A is the normalization factor), we obtain 
that E(k)dk = § ^2A 3l2 dt k . Now, A = \ 3 from the normal- 
ization of the turbulent spectrum (Eq. k^ = (ViA) -1 with 
Vl the velocity of the largest eddy, and V\ = |V? also by nor- 
malizing the power spectrum (see CH03). We then end up with 



1 V 2 

E(k)dk = -—dt k , 

2 ?L 



(11) 



which can be substituted into all the integrals, putting them into 
a simpler form. For instance, Eq. (O becomes for particle i 



V: 



«l Jt, { \ti + h) 

Similarly, the cross term becomes 

. vl 



2\ 



tL h 



2U f' L 
' + f 2 Jru 



dtk 



1 



ti 



V 2 

g 



2\ 



tr 



tk + 



ti + t k 



ti + t k 



2t t 



t L h + 1 2 



tk + 



ti + t k 



(12) 



(13) 



(14) 



With t* 2 = max(f*,fp and f, ; < t* 2 < ?l since t* refers to an 
eddy's turn-over time. We now solve for AV 2 2 by splitting the 
integral in Eq. ( TTZl ) at t* 2 and subtracting the corresponding V c , 
terms from Eq. ( fT~4-b to get 



V 2 

AV 2 2 = -i 

'L 



tk + 



h + t k 
2h 



tk + 



h+t k 



t\ +t 2 





tL \ 


tk + — l — 


+ (1^2) 


h+tk 





, (15) 



V 1 

AV 2 2 = -i 



t"\ 



tk + 



h+t k 



h-h 
h +t 2 



' i 



h +t k 



+ (1 <-»2) 
(16) 



= AVl + AV 2 



This is perhaps the most concise way to write the expressions for 
AV 2 2 . The first term we call AVj] since this term involves class II 
(fast) eddies. If t* 2 = ?l (heavy particles) all eddies are fast and 
only this term remains. Conversely, if t* 2 = t^ (small particles) 
the contribution from AVn vanishes and the second term, AVi, 
determines relative velocities. In the intermediate regime, t n < 
fj 2 < ?l, both terms contribute. Written in terms of the Stokes 
numbers these terms becomes 



AVf = 



h-h 



= V 



tL h + h 
2 Sti-St 2 



h +t k 



E St! + St 2 



Stt 



+ (1^2) 



Stf 



St* 2 + Sti 1 + St! 



-(1^2) 



(17) 



V, 



AVs = — 



tk + 



l 



h +h 



+ (1^2) 



Vi 



(St: 



Re- 1 ' 2 ) + 



St 



—rr + (1 2) 

St!+St* 2 Sti + Re- 1 ' 2 ') 



(18) 



Note again that since t n < t* 2 < f L we also have that Re~ 1 ^ 2 < 
Stj 2 < 1. Below, we will first solve for Stj 2 , and then consider 
solutions for AVn in various limiting cases of the particle stop- 
ping times. 

3.2.1. Solving for f 

The relative velocity b etween a particle with stopping time f s and 
an eddy k, is given bv lVolk et all (1 1980b . Eq. (15): 



Vrdk) 2 = V 2 +2 f dk'E(k')(-^ 

Jk L \t s + 



(19) 



V D is any systematic velocity component not driven by turbu- 
lence - such as due to pressure-gradient driven azimuthal head- 
wind, the ensuing radial drift, or vertical settling under solar 
gravity. We can integrate this equation in the same fashion as 
Eq. ( fT4b and arrive at 



V 2 

v Ie m 2 = v 2 + ^ 

tL 



t: 



h + tk 



= Vl + 



tL 



1 



1 



1 +VL) 

(20) 



in which y - t k lt s . Also, using the definition for t k (see text 
above Eq. (fTTTi). k* can be expressed as (k*) 2 = (2A)~ 3/2 t*~ 3 = 
|y~ 2 t\f ~ 3 . Inserting the expressions for k* and V 2 A into Eq. (O, 
assuming that V — for simplicity (see however Sect. 13.31 ). we 
obtain: 



?r? l=,; ' 



(21a) 



4 



C. W. Ormel and J. N. Cuzzi: Closed-form expressions for turbulent relative velocities (RN) 




Fig. 1. The function \y*{y* - l) 2 - 1/(1 + y*). If f s <k f L and 
V — (no systematic velocity drifts; see Sect. |3.3l ) this equation 
is equal to zero and we find a solution y* = f/t s ~ yl - 1.6. On 
the other hand, for f s ~ f L , the RHS of Eq. ( |21d) is « -0.5 and 

y* * 1. 




2 * / * 
3# 0/ 



l) 2 



1 + */* 1 + !/L ' 



(21d) 



where we have defined i/* = f*/f s and j/l = t^/t s = St . The 
LHS of Eq. ( I21dl ) is pl otted in Fig. [T] If j/ L » 1 for small par- 
ticles, the RHS of Eq. (12 1 db is negligible and the numerical so- 
lution for y* becomes y* ~ y* a - 1.6, or f* w 1.6f s . On the other 
hand, when f s nears ?l, the -1/(1 + yi) term causes the RHS 
of Eq. ( I21dl > to drop to -0.5, and y* — > 1. For f s > ?l we al- 
ways have that f* = ?l; '•<?., for such a particle all eddies are of 
class 2. In Fig.|2]we compare the exact solution (dashed line) for 
f with the f « i/*? s = 1 .6f s approximation (solid line; in both 
cases t* < ?l is simply enforced), and the empirical function 
k*/k L = 1 + ±Sr 3/2 (dotted line; see Sect. ED), 

The exact solution for AVn (Eq.fToT) is given in Fig.[3]both 
for fi » ?2 (solid curve) and for particles of equal stopping 
times (dashed curve). A Reynolds number of Re = 10 8 has been 
adopted. 

3.3. The role of V a : eddy-crossing effects 

Systematic velocities V due to vertical settling , and pressure- 
gradie nt headwinds and drifts, will occur (eg. Nakagawa et al.l 
( 1986)). Because particles drift through eddies, their transit time 
is affected (because V re i is larger) and the bo undary between 
class I and II eddies shifts. Cuz zi et al. I (11993b include this ef- 
fect, due to vertical settling, in their model of particle diffusion 
(their Eq. (43)). The model presented here offers a generalized 
way of treating this effect, which we will only sketch here. 

Repeating the procedure outlined in Sect. 13.2. 1 I but retaining 
the Vo term in V re i (Eq. (l20l). we end up with Eq. ( I21db including 
a correction term 



y (y - 1) 



l +y* 



= F(y) = - 



St 



l y ( 2 

1 +st + si ~v}, 



-4, (22) 



where we have substituted St = 1/i/L- The correction term 
can be roughly constrained using an estimate of the system- 
atic drift velocity V a ~ (St/(St + l))/3V K , where V K is the 
Keplerian velocity at distance a from the Sun, D. is the orbit fre- 
quency, and p = {H g /a) 2 is a radial pressure gradient parameter; 
also we take V B = a 1/2 c e (see, eg.. iNakagawa et ail d 19861) or 



ICuzzi & Weide nschilling (2006b). Then 



Vo = st pv K 

V g St + 1 a l ' 2 c g 
St 13 



St + 1 a 1 ' 2 ^' 2 
and Eq. d22l becomes, 



St /3aQ 
St + 1 a l l 2 H g Q. 

- St (£\ 1/2 

~ St+ 1 \a) ' 



(23) 



Fig. 2. Three different assumptions for f (or the related k*) are 
shown here. 



f^-L] =(2.1) -if 



1 



rel 



f (y* - l) 2 = \U . . 

- \ 1 + y* l+y L 



1 



(21b) 



(21c) 



F(y) = 



St / p/a 



1 + St \ 1 + St 



- 1 



(24) 



Normally P ~ 2 x 10 -3 is assumed dNakagawa et all [l986; 
ICuzzi et alll 1993b . but its real value, and that of a, are not well 
known. Equation (l24l shows that for a given value of St, F(y*) 
increases with increasing ft /a. Consequently, y* = f/t s is also 
higher (see Fig. [1). The boundary between the class I and II 
eddies therefore shifts to higher values of f, that is, there are 
less class I eddies for high pi a and the St* = 1 upper limit 
(when f = t]_) is reached at lower Stokes numbers. Inserting 
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i i iiiiiii — i 1 1 mill — i i iiiiiii — i 1 1 nun — i 1 1 mm i i in 




- - h=t 2 
h > h 

Fit Eqs. (27, 28, 29) 



I I I I Mill I I I I 



thus 



1(T 5 1(T 4 1(T 3 1(T 2 lfT 1 10° 10 1 
Sti =h/t L 

Fig. 3. Exact solution to Eq. dT6b for AVn in the case of identical 
particles (dashed line) and t\ » f 2 (solid line) for a Reynold 
number of 10 8 . The dotted curves are approximations to Eq. ([TrjT i 
given by Eqs. 



the definition of F(y*) (LHS of Eq. d22>) into Eq. (|24| with 
y* = t^fts = St~' and solving for St, we find that the Stokes 
number at which St* = 1 occurs at 



St s ,.=i = 




(25) 



For example, for f}/a = 1, St* reaches its upper limit at St « 
0.67. 

In the small particle regime (St <K 1), however, the ex- 
act value of pi a is unimportant since F(y*) is always close 
to zero, and the j/* approximation is justified. It is only for 
f3/a > Sr 1 that the RHS of Eq. d24]i starts to becomes signif- 
icant and {/*>{/*• This is the weakly-turbulent or non-turbulent 
regime where class II eddies dominate even for small particles. 
In practise, however, it means that eddy crossing effects are im- 
portant only if turbulence is very weak and we will not treat them 
further in this paper. 

3.4. Limiting solutions 

As intuition-building examples we obtain simple, closed-form 
expressions for AV 2 2 in various limiting regimes from the f-space 
solutions; similar results are easily obtained from the fc-space 
solutions (Sect. |3TTT >. Without loss of generality we take particle 
1 to have the largest stopping time, i.e., t\ > f 2 and f* 2 = f*. 
Moreover, we assume that t n <K t^; i.e., that Re 1 ^ 2 » 1 and 
there is an extended inertial range of eddies. Recall again that 
St* 2 = Re~ 112 for t\ < t 1] jy\, and that St* 2 will not exceed 1. 



3.4.1 . Tightly coupled particles, h,h< t n 

In this limit all eddies are of class I and AV 2 2 
particle, the second term on the RHS of Eq. ( fT7] l is negligible 



AV 2 . For each 



AV: 



12 



V: 



2 Sti -St 2 
; Stj + St 2 



St 



St* 



St, +Re- 1 ' 2 Sb+Re- 1 ' 2 



(26) 



In the very small particle regime (fi «: f 7 ), St, <s t n /tL = Re 1/2 
and 



av 2 2 = y*-(Sti 



St 2 ) 2 



(27) 



Since V 2 - |V 2 /?e^ 2 - jV^t^/t^, this expression transforms 

directly to AV; 2 = V3/2(fi - f 2 )V^/f^, in good agreement with 
the heuristic, al though physically motiv ated, expression AVi 2 = 
V,,(fi - t 2 )/t v of lWeidenschillin j dl984 . 

3.4.2. Intermediate regime, t n < t\ < f L . 

If t\ (the stopping time of the larger particle) approaches the 
Kolmogorov scale, two changes occur. First, the Stj/(Stj 2 + Sti) 
term in Eq. (jT7j now becomes linear with Sti, since St* 2 grows 
proportional to St; (the second term is still negligible throughout 
most of this regime). Relative velocities therefore increase as the 
square-root of stopping time. Second, class II eddies also con- 
tribute to AV 2 2 (Eq. [T8l>. This contribution scales also with Sti, 
but is significantly larger and does not disappear when t\ — f 2 . 
From a physical point of view, class II eddies act as small, ran- 
dom kicks to the particle trajectory, while two particles captured 
by a class I eddy are subject to the same, systematic, change 
in motion. Class II eddies are therefore much more effective in 
generating velocity differences for similar-sized particles. 

In the "fully intermediate regime", i.e., t n « f| < ?l, we 
can also ignore the Re~ 1 ^ 2 terms in Eq. ([TBI . In addition, the 
t* lh = yl approximation holds. Upon writing St 2 = eSti, Eqs. 
( fT71[T8l l become linear with Sti and we can write AV 2 2 as (see 
Appendix iBt 



AVf, 



Vi 



2j/ a - (1 + e) + 



1 



l+e\l+{/a y a + e/ 



Sti, (28) 



where e < 1 is the ratio between the stopping times and 
i/ a = 1.6. For fi » t 2 we then find that AV 2 2 * 3.0V 2 Sti, 
while for equal particles the numerical factor goes down to 2.0. 
Written in terms of stopping times the relative velocities be- 
come , AVi 2 = [1.7 -f 2,1"1V l V f i/ f L- This also compares well 
wi thlMidenschilling (1984) fits for this regime (who gives pre- 
factors of 2.1 and 3.0, respectively). Note, however, that our full 
expressions for AV (Eqs. [16] [17] [T8l also capture the behavior 
near the t n and "turning points" (see Fig. [3). 



3.4.3. Heavy particles, fi > t h 

If fi > ?l, St* 2 = 1 and there is no contribution from class I 
eddies (Eq.fTTIi. Also, we can neglect the Re~ 1 ^ 2 terms in Eq. ( fT8l 
and the relative velocities simply become 



AV 2 = AV 2 = Vi 



1 



1 



1 + Sti 1 + St 2 



(29) 



This result can, of course, directly be obtained from the 
V P i terms (Eq. ( fT2l ) since the cross-term vanishes in this 
regime. For small St 2 relative velocities are still ~ V g ; 
however, if both Stokes numbers are large, the relative 
velocity decreases roughly with the square root of the 
smallest particle stopping time. Note that the linear fit of 
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< (A n 6v) 2 >yj, / < Sv 2 > 




log St, 



Fig. 4. Contour plots of particle-particle, turbulence induced, relative velocities AVn normalized to V g . -A- Numerical results of 
Mar kiewicz et al.l (Q991), without inner scale (Re — > oo). -B- Analogous result from our closed-form expressions with the fixed 
y* » y* = 1.6 approximation (Sect. !3.2.Tb . -C- Like B, but with an exact solution for y* and with Re = 10 s . -D- Using the CH03 
formula for k*, k* jk^ = 0.5St* + 1, and also with Re = 10 8 . Contours are drawn twice per logarithmic decade (at AVn/Vg = 3x10' 
and at 10') with an additional contour at 0.8 and 1.15. 



IWeidenschiiiingl d 1984 in this regime (his Eq. (15)) is inappro - 
priate (see, however, IVolk et alJfl980h:IWeidenschijlingl ([l988); 
IWeidenschilling & CuzzH d!993l) : ICuzzi et all d 1993b in which a 
square-root fal l off is advocated ). Since an explicit, closed-form 
solution to the lVolk et al](fl980l) and MMV expressions for AV12 
has n ot previously been availab l e, many dust coagulation model s 
(e.g.. ISuttner & Yorkd d200ll) ; iDullemond & Dominikl J~2005): 
lOrmel et all d2007l) ) have rehed on the IWeidenschiiiingl 7l984) 
fits to calculate relative velocities. Turbulent motions and rela- 
tive velocities for particles in the f s > ?l regime have therefore 
been underestimated in these calculations. However, concerning 
these works, we also think no major conclusions have been af- 
fected, since the error is introduced only for large dust particles, 
that is, if the system is already well evolved. 

3.5. Contour plots 



Following I Volk etldl (Q~980) and MMV we also present our re- 
sults as contour plots. Figure HJ\ shows, for comparison, the 
results of MMV, obtained by numerical evaluation of the inte- 



grals involved without an inner turbulent scale [Re — > 00). The 
next three panels of Fig. [4] show the result using our closed-form 
expressions derived from Eq. dToT ). In panel B, the y* approxi- 
mation has been used and, like Fig. 2 of MMV (panel A), the 
inner scale of the turbulence is extended to infinity so that Eqs. 
d28ll29l ) apply. Somewhat systematically higher values for AVn 
when compared to MMV can be explained by the CH03 approx- 
imation for V p (see Eq. (|6]l) but these discrepancies are less than 
~10%. In panels C and D we show the contour plots correspond- 
ing to the other formulations for k* (see Fig. |2), i.e., the exact 
solution for y* (panel C) and the CH03 empirical approxima- 
tion (panel D). The differences between these three methods for 
determining k* differ around the St = 1 point (see Fig. |2) and 
are reflected in the contour plots. For St » 1, panel C compares 
best to the numerical result of MMV, but no significant errors are 
made when using the y* approximation or the CH03 formula for 
k*. 

In panels C and D of Fig. |H a Reynolds number of Re — 10 8 
has been adopted. For St < 10~ 4 , therefore, velocities are greatly 
suppressed since only class I eddies remain to generate relative 
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velocities and relative velocities disappear completely for equal 
friction times. Also, the contours are much closer spaced since 
in this limit the velocity AVn is proportional to St (see Eq. (1261)). 



4. Conclusions 

We have exten d ed and , essentially, completed the work of 
ICuzzi & Hoganl (|2003), who derived explicit, closed-form 
expressions for particle velocities i n turbulenc e base d on 
the physics orig i nally developed by IVolk et aTl dl980t) and 
Mar kiewicz et al.l £1991). Within the framework of this physics, 
the only approximations used here are in Eq. (O for the particle 
velocities (where a posteriori comparisons with exact numeri- 
cal solutions indicate the approximation is well justified) and in 
Eqs. d20j > et seq where the systematic velocity V is neglected 
to simplify calculating the boundary between eddy classes (gen- 
eralizing this step should be straightforward, however). The full 
analytic expression for AVn is given by Eq. ( fT6] l (or by the sum 
of Eqs. ( fT7l[T8] l), but more simple, explicit expressions apply in 
restricted regimes (provided Re 1/2 » 1): 

- Equation d27| >, in the very small particle limit {t\ <& t n ); 

- Equation d28l ). in the "fully intermediate" regime, i.e., for 
t n <k h «: f L ; 

- Equation j29t . for t\ > t^. 

Near the t\ = t n and t\ = ?l turning points the behavior is 
more complex (see Fig. [3]) and for accurate analytical approx- 
imations one has to revert to the full expressions for AV given 
by Eqs. CSlIIllIIl. 



10 u = 



10"' = 



mill i i ii 




10° 10 1 10 2 10 3 10 4 10 5 10 6 

k/k L 

Fig. A.l. The functions g (solid lines) and h (dashed lines) plot- 
ted for three different Stokes numbers at a Reynolds number of 
10 8 . After k > k* the functions show power-law behavior. The 
power spectrum (weighted by k) is plotted by the dotted line. 
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Appendix A: A more accurate closed-form solution 
for V p and all related velocities, using power-law 
approximations to the functions g and h 

In Sect. |2 the very simple approximation g(x) = h(x) = 1 was 
introduced for all stokes numbers St and eddy scales k. It proves 
to be quite adequate for most purposes; however, as noted in 
Sect. 13.51 small inaccuracies remain at the 10% level because 
the approximation overestimates the contributions of fast eddies 



to Vp and other velocity components. Figure IaTI shows the de- 
tailed behavior of the functions g and h for Re — 10 s for Stokes 
numbers of St = 0.01, 0.1, and 1.0. The inflection point for all 
three values of St is at k = k* (recall that k*/k L * 1 + ±Sr 3/2 ). 
For k > k*, the functions are well approximated by power- 
laws of -1/3 and -3/4, respectively, i.e., g(k) = (k/k*y 1 ^ 3 and 
h(k) - (k/k*y 3 ^ 4 . The success of the approximation of Sect. [2] 
is due to the fact that the power in the weighting function kE(k) 
(dotted line; we multiply with k since we compare logarithmi- 
cally) decreases rapidly with increasing k; thus by the time the 
assumption g(k) = h(k) = 1 becomes really bad, the relative con- 
tribution of successive terms has become small. For small St, the 
weighted contribution of eddy power has already become very 
small even before k ~ k* (the logic of CH03). For St = 1 or 
larger, the weighting function has dropped by nearly an order of 
magnitude by the time h(k) (the faster-decreasing function) has 
dropped to 0.3 (dashed-dotted line), and this seems to account 
for the success of our simple assumption. 

This behavior can be understood from the definition of x — 
KtkkV K \. For k » k*, K as 1 and V re i < V g are both constant. 
Then, because oc k~ 2/i , x scales as x 00 k x ^ and becomes 
large at large k. Since g(x) = arctan(y)/^ oc x~ l for large x, 
we get that g(k) oc k~^ 3 . Similarly, h(x) oc AT 2 / 3 , which is a bit 
shallower than the -3/4 exponent observed over most regions of 
interest (Fig lA.lt . While the -2/3 exponent is reached at large 
k, the -3/4 exponent seems more appropriate at intermediate k. 
Yet, in our subsequent analysis, we will use the large-£ limit for 
this exponent (-2/3) because it simplifies the math. Thus, we 
approximate the g and h behavior as follows: unity for k < k*, 
and power laws in k/k* with exponents of -1/3 and -2/3 for 
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k > k*. Then Eq. © becomes 
V p 2 = dk2E(k)(l-K 2 ) 

Xk„ 
dk2E(k){\ -K) 



kY U3 1^ 
k*J +K \k* 



(A.l) 



Where we have that k n < k* < k^ such that in the case of very 
small or very large particles one of the integrals vanishes (Sect. 
13.2. U . Since the approximation g = h = 1 still holds for k < k* 
(or for t > f) the velocities resulting due to class 1 eddies (Eq. 
fFft ) are not affected; the new approximation only affects Eq. 
<H1J. By writing K = St/(St + ;T 2/3 ), E(k) = E L (k/k L y 5/3 oc 
jT 5 / 3 with x = k/k^ the solution to Eq. jA.lj involves integrals 
of the form 



dxx- 5/3+l) 



St 



\St + x- 2 / 3 



n=[l,2] 



with n = 2 for the K 2 term and p = -1/3 or -2/3. These in- 
tegrals can be solved analytically. Going to "t-space", however, 
gives somewhat cleaner solutions and we will from here on fol- 
low that approach and show how it affects relative velocities, 
i.e., AV n - After the change of variables faff = (k/k*)~ 2/3 ) the 
second term of Eq. ( 1A.U becomes 

ir^ (i -^(?r +(i -^(?)- (Aj) 



We now introduce the dimensionless variable y = t^/t s (cf. Eq. 
(120b). Then t k /f = y/y* with y* = f/t s . Also K = 1/(1 + y) and 
1 - K = y/(l + y) and Eq. ( IA.3b becomes 



(A.4) 



V t r> is „-i/2 ,,2 

» s \ a i n-i/2 f , / n-i f 
— dy{y) 1 —— + (y ) 2 

= ^ {(y*r 1/2 [m)f t ^ + (y*r l [m] f t ^} (a.5) 

in which the functions h(y) and I g (y) are defined as 

r» z 3/2 [2 \ 

hiy) = J Q dz ^ = I - 2 1 y/y + 2 arctan( V^) (A.6) 

l g {y)= [ J dz-^— 2 =^±l^-21og(l+y) (A.7) 
Jo (1+z) 2 1 + 1/ 



This gives us several terms that contribute to AVn. Collecting 
these terms, the new expression for AVn becomes 



-1/2 



+ " 7^ 



4(z/) 



t\lh 



t n lh 
t"Jh 



+ (1 <->2) 



+ (1«2) 



(fi > fe). 



(A.8) 



Although still fully analytical, this more accurate expression for 
AVn is also more complicated and we did not present it in the 
main body of the paper. Equation ( 1A.8I ) is useful, however, for 
readers whose applications demand this higher level of accuracy. 

Appendix B: Derivation of Eq. d28b 

(A.2) 

We consider the limiting case of t n <K t\ «: ?l. The y* approxi- 



mation for Stf 2 then holds, i.e.. 



St* 2l « i/*Sti with y a = 1.6. We 



will now argue that we can neglect the Re~ l t 2 terms in Eq. (IT8l . 
For particle 1 this is obvious since Sti » Re~ l t 2 . The last term 
(where Re^ 1 ^ 2 is in the denominator) then becomes simply -Sti . 
However, for the interchange term a similar approximation 



St 2 + Re 



-1/2 



St 2 , 



(B.l) 



is not that obvious since we have not put a constraint on St2. 
For example, if St2 <K Re~ 1 ^ 2 the Re ' 2 term dominates the 
denominator. However, in that case this term and its approxi- 
mation are small anyway compared to —Sti, such that by mak- 
ing the approximation in Eq. flB.lt our final result is not af- 
fected. Similarly, if St2 ~ Re~ 1 ^ 2 , Eq. (IB. Il l (which goes to 
~ jRe~ l l 2 ) or its approximation (~ Re~ 1 / 2 ) are insignificant 
since Sti » Re' 112 . Only if e ~ 1, i.e., St 2 » Re' 11 2 , do es 
the St2 term matter, but then the approximation in Eq. flB.lt is 
well justified. All terms in Eq. (IT~8b are then linear in Stokes and 
we can reduce it to 



A ^n / 1 e 2 



1 + yt e + yl 



St, 



with e = St 2 /Sti < 1. Similarly, Eq. dTTb becomes 



AVf 1 - e / 1 



l+e\y* a + l y* + e 



Sti. 



(B.2) 



(B.3) 



Combining these expressions and collecting the 1/(1 + y* a ) and 
e 2 /(y* a + e) terms then gives Eq. 



The expressions for AVn now consist of several contributions. 
First, Ihiy) and I g (y) are evaluated at both the upper (y*) and 
lower (y n ) limits. This must be done for both particles 1 and 2, 
because the AVn term (Eq. ( fT8l ) has separate contributions from 
each particle. For the particle of highest friction time (say this 
is t\) the power-law approximation for g and h holds over the 
range AVn is calculated, i.e., t s < fi* < t\ 2 = t*. However, for 
the second particle the power-law approximation only holds for 
hk < t* v while for the remaining range over which the integral in 
AVn is evaluated, i.e., f 2 < tn < f*, the g — h — 1 approximation 
applies. 



